Excitation of multiple giant dipole resonances: 
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Abstract 

The effect of deformation on the excitation of multiple giant dipole resonances is 
studied. Analytical expressions are derived in the framework of the interacting bo- 
son model for the energies and El properties of giant dipole resonances in spherical 



O I and deformed nuclei, and a numerical treatment of transitional nuclei is proposed. 
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■ ' Coulomb-excitation cross sections are calculated in ^38^ g^^^^^ ^^ ^^le samarium iso- 



topes. 



Nuclear multiphonons both at low and at high energies have attracted much 
interest lately [1-4]. At low energy, the controversy centers around the collec- 
tive character of multiphonon states. While in vibrational nuclei firm experi- 
mental evidence now exists for states with triple quadrupole phonon character 
(see, e.g., [5]), in deformed nuclei the collective character of double f3 and 7 
vibrations is still a matter of acrimonious debate. At high energy, the study of 
multiphonons has pointed out the limitations of the small-amplitude approx- 
imation for vibrational collective motion and of the linear approximation for 
the external exciting field. These assumptions are routinely made for the first 
phonon but recent studies of double phonons [6-8] have shown the importance 
of anharmonicities and non-linearities in the excitation of large-amplitude vi- 
brations. Similarly, anharmonicities and non-linearities in nuclear vibrations 
at low energy have been shown to play a crucial role in the calculation of the 
heavy- ion fusion cross sections at energies close to the Coulomb barrier [9,10]. 
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In this letter yet another aspect of phonon excitations at high energy is investi- 
gated, namely the modification of the excitation of the giant dipole resonance 
(GDR) away from shell closure as a result of deformation. Already at the 
level of the single GDR the strength function shows a splitting into two peaks 
associated with the two different frequencies of vibration along with or perpen- 
dicular to the axis of axial symmetry. The question under scrutiny here is how 
deformation infiuences the El strength to the double GDR (DGDR) and what 
its effect is on the Coulomb excitation cross section. The proposed approach 
makes use of the interacting boson model (IBM) [11] for the description of the 
low-energy collective levels. These are coupled to the GDR excitations mod- 
elled as p bosons. The advantage of using the IBM is that the multiphonon 
states are obtained as exact eigensolutions of the hamiltonian and, therefore, 
no folding procedure is required to obtain the double phonon states, as is 
for instance the case in [12,13]. The folding procedure is only approximately 
correct for vibrational or well-deformed nuclei (in the latter case the folding 
must be done in the intrinsic frame assuming additivity of intrinsic phonons) 
but there is no simple recipe in the intermediate case of transitional nuclei. 
Another advantage of the IBM is that calculations are quick and this enables 
an easy estimate of the excitation cross section to the triple GDR concerning 
which experiments are currently planned [14]. 

A general and appropriate basis to discuss the problem of multiple GDR ex- 
citations is of the form 

\aLLxnR;JMj). (1) 



The multiple GDR is built on a low-energy nuclear state characterised by 
a^L where L is the angular momentum of the state and ai any other label. 
The multiplicity of the GDR is indicated by n (i.e., n = 1 for a single GDR, 
n = 2 for a double GDR, etc.) and its angular momentum by R. The single 
resonance is approximated as a p boson; the allowed angular momenta of the 
multiple GDR are i? = n, n — 2, . . . , or 1. The basis (1) can be referred to 
as a weak-coupling basis in the sense that the angular momenta L and R are 
good quantum numbers and are coupled to total angular momentum J with 
projection Mj. The weak-coupling basis arises naturally when the interaction 
between the GDR and the low-lying states is weak or is of dipole type L ■ R. 
More general nuclear interactions are not necessarily diagonal in the basis (1). 
Specifically, the interaction predominantly responsible for the splitting of the 
GDR in deformed nuclei is of quadrupole type Ql ■ Qr and is not diagonal in 
the weak-coupling basis (1). From the analogous problem in the particle-core 
coupling model [15] it is known that the diagonalisation of Ql ■ Qr in the 
basis (1) gives rise to a strong- coupling basis with the quantum numbers [16] 

\aLXn;ajKjJMj). (2) 



The angular momenta L and R no longer are conserved quantities and are 
replaced by Kj, the projection of the total angular momentum J on the axis 
of axial symmetry. The choice of basis depends on the competition between 
various terms in the nuclear hamiltonian. Large splittings in L or i? induce 
the weak-coupling basis (1); if, in contrast, these are small in comparison with 
the quadrupole coupling between the low-energy states and the GDR, such as 
is the case in well-deformed nuclei, the strong-coupling basis (2) is obtained. 

The above remarks are rather general and model independent. Low-energy 
nuclear states and their coupling with the GDR can be modelled in several 
different ways and a convenient one is in the context of the IBM [11,17]. In 
deformed nuclei the problem was worked out analytically for a single GDR by 
Rowe and lachello [18], while numerical results were presented for the single 
GDR in [19-22] and later for the double GDR in [23] . In this letter analytical 
results for the energies and El transitions in spherical and deformed nuclei 
are generalised to the multiple GDR, and numerical results are presented for 
intermediate cases. From these results the Coulomb-excitation cross sections 
are calculated taking into account the dynamics of the collision, which was 
lacking in previous IBM treatments of the single and the double GDR. 

It is assumed in the IBM that collective nuclear states can be described in 
terms of A^ s and d bosons where A^ is half the number of valence nucleons [11]. 
The dynamical algebra of the model is U(6) in the sense that a single of its 
representations (namely the symmetric one, [A^]) is assumed to contain all low- 
energy collective nuclear states. To this space are coupled the multiple GDR 
excitations. Assuming that a single GDR is described by a p boson, multiple 
GDR excitations are represented by the direct sum (0, 0)©(1, 0)©(2, 0)©- ■ ■ of 
symmetric representations (n, 0) of U(3). The dynamical algebra of the coupled 
system is thus U(6) © U(3) with the proviso that several U(3) representations 
must be taken to build the model space. 

The model hamiltonian has the generic form [17] 

H = Hsd + Hp + Vsd-p- (3) 



The first term is the usual IBM hamiltonian [11] which gives an adequate 
description of the low-energy spectrum of spherical, deformed and transitional 
nuclei. The second term in (3) governs the multiple GDR spectrum and is of 
the form 

Hp = eph + aph{h + 3) + PpR^ , (4) 



where n is the p-boson number operator and R the associated angular mo- 
mentum operator. The coefficient e^ represents the unperturbed single GDR 



energy. The second term in (4) induces a diagonal anharmonicity in the excita- 
tion energy of the multiple GDR. The interaction term in (3), finally, acquires 
the form 

Vsd~p = aofidfi + 2aiL ■ R + 2V3a2QL ■ Qr, (5) 



with a monopole, dipole and quadrupole interaction. More complicated, higher- 
order interactions that mix excitations with different n are not included here. 
The hamiltonian (3) can be diagonalised numerically in the U(6)®U(3) model 
space. For particular values of the parameters analytical solutions are avail- 
able. Two classes of analytically solvable hamiltonians exist. 

(i) Weak-coupling limit. A weak-coupling hamiltonian can be defined for each 
of the three limits [U(5), SU(3) and 0(6)] of the IBM and is obtained for 
«o = 0^2 = in (5). In the particular case of the U(5) limit, only a2 = is 
required. In the U(5) weak-coupling limit states are labelled by 

I [N]ndVnAL x (n, 0)R; JMj) , (6) 



where rid is the number of d bosons, v the rf-boson seniority and n^ an addi- 
tional quantum number related to the pairing of triplets of d bosons [11,25]. 
The energy eigenvalues of the states (6) are 

E = edUd + adUdiud + A) + 'ydv{v + 3) + {(3d - ai)L{L + 1) + 

epU + apn{n + 3) + {(3^ - ai)R{R + 1) + aQUdU + aiJ{J + 1), (7) 

where 13d is the parameter associated with U' in Hgd- 

(a) Strong- coupling rotational limit. The appropriate classification of the low- 
energy states in this case is SU(3). This symmetry requires that the quadrupole 
operator Ql be a generator of SU(3) and that ao = and Pp = Pd = cti — f ^2, 
In the SU(3) strong-coupling limit states are labelled by 

\[N]{\sd, fisd) X (n,0);(A,/i)KjJMj). (8) 



The labels {Xsd, l^-sd) are associated with the SU(3) algebra of the s and d 
bosons. They characterise the band structure of the low-energy spectrum; for 
example, the ground band has {Xsd, fJ-sd) = (2A^, 0). The energy eigenvalues of 
the states (8) are 



E= {asd - a2)[Kd{Kd + 3) + ^isd{fJ-sd + 3) + XsdfJ'sd] + 
epH + {ap - a2)n{n + 3) + 
a2[A(A + 3) + /i(^ + 3) + Xfj] + (en - fas) ^(^ + !)• (9) 



The splitting of the GDR comes about because of the fourth term in (9). The 
(GDR)" excitation sphts into n + 1 peaks corresponding to (A, fi) values 

(A, ^) = (2iV + n, 0), (2A^ + n - 2, 1), . . . , {2N -n,n), (10) 



where 2N > n is assumed. The energies of the different peaks are found 
from (9). The values of Kj and J allowed for a given (A, fi) representation are 
given by Elliott's rule [11,24]. 

The Coulomb excitation of GDRs occurs predominantly through El. In the 
context of the present model an El excitation corresponds to the creation 
of a p boson (annihilation in case of El de-excitation) and thus the electric 
multipole operator A4(Elfi) [15] is parametrised as C{p^^+p^i)- The calculation 
of El transition probabilities requires the matrix elements oip^ in the basis (6) 
or (8) which can be done by standard group-theoretical techniques. Analytical 
expressions are found in the two limiting cases. 

(i) Weak- coupling limit. For the GDR excitations built on the 0+ ground state 
results up to the DGDR are shown in Fig. la. Generally, for the B{E1) values 
between multipole GDR excitations built on the 0"*" ground state one recovers 
the independent-quanta result [15] 



J2 B{E1; nd = Ox (n, 0)Ri = J,; J, ^ na = x {n - 1, 0)Rf = Jf, Jf) 
f 

= ^ B{E1; anJi -^ a{n - 1) J/) = nB{El; an = l^an = 0), (11) 

/ 

where a denotes all other quantum numbers that cannot change. 

(a) Strong- coupling rotational limit. Results up to the DGDR are shown in 
Fig. lb. Although all individual i?(El)s are known, for simplicity of presen- 
tation only the summed strengths J2k'L' B{E1; {X,fi)KL -^ {X', ii')K'L') are 
shown. In the limit of large boson number A^ one recovers harmonic results 
that have a simple geometric interpretation. For example, the B{E1) values 
from the 0^ ground state to the two 1^ GDRs are 1 and 2, respectively, the first 
associated with an oscillation along the axis of symmetry (say the z direction) 
and the second with oscillations in the x and y directions. For the single-to- 
double GDR excitation one finds B{E1) values which are, for N -^ oo, 2, 2, 1 
and 3. The large-A results can be generalised to (GDR)". 

Coulomb excitation in heavy-ion collisions is usually described by treating the 
relative motion classically while the internal degrees of freedom of the colliding 
nuclei are accounted for quantum mechanically. The operator responsible for 
the excitation depends on time through the relative distance. For relativistic 
collisions its expression is as in equation (35) of [7] where each term of the 



multiple expansion of the external field factorizes into two elements. The first 
depends on the collision properties, the second on the structure of the nucleus 
being excited. In the present study only contributions from the A^(El) matrix 
elements are considered and those are calculated within the model described 
above. The solution of the time-dependent Schrodinger equation leads to a 
set of coupled differential equations for the probability amplitudes to excite 
the (GDR)" states. For each impact parameter b these equations are inte- 
grated along the appropriate classical trajectory. For each (GDR)" state the 
total inelastic cross section is then obtained by integrating the correspond- 
ing probability over all impact parameters, starting from a mininum value 
Vin = IMIA^^ + aI^^ - 0.75{A^^^^ + .4^^/^)] fm [26]. 

The above formalism will now be applied to the relativistic Coulomb excitation 
of the single and double GDRs in ^^^U and in the chain of even isotopes ^^^Sm 
to ^^^Sm. The former is an example of a well-deformed nucleus while the samar- 
ium isotopes exhibit a change from vibrational to deformed as A increases. The 
description of such structural changes requires the use of a transitional IBM 
hamiltonian and hence the following analysis is not confined to any of the 
previously discussed analytical limits but always involves a numerical diago- 
nalisation. For ^^^U, a consistent-Q [27] hamiltonian Hgd = kQ^-Q^ + k L-L is 
used with k = —16 keV, k' = 1.5 keV and x = —0.72. These parameters yield 
an adequate description of the ground-gamma band splitting, of the moments 
of inertia and of the E2 transitions from gamma to ground band. The Hgd 
hamiltonian for the Sm isotopes is taken as in [21]. The additional parameters 
ep, ttO) ^2 and ( in Hp, Vsd~p and in the El operator are given in Table 1. 
They have been chosen as to reproduce the observed photoabsorption cross 
section [28,29] to the first GDR. Agreement is obtained if to each eigenstate 
is associated a spreading width Fj = 0.007Ef-^ in ^^^U and Fj = 0.029£'/-^^ in 
the Sm isotopes (with Ei and Fj in MeV). To reproduce the photoabsorption 
cross section of a well-deformed nucleus one uses the fact that the energy split- 
ting of the GDR is very sensitive to the parameter 02- Then, once 0:2 is fixed, 
one varies ao to sligthly change the contribution of each peak to the energy- 
weigthed sum rule and tp to shift the energy of the dipole states. Finally, the 
parameter ( is obtained from a global normalisation. 

Recently, an experiment was done for ^ssy _|_ 208p|^ ^^ g^g GeV/A (the data 
analysis is in progress [30]). The result of the corresponding calculation is 
shown in Fig. 2. The full line corresponds to the total cross section obtained 
by smoothing the cross section to each discrete state with a Lorentzian having 
a width of Fi = 2.5 MeV and F2 = v^Fi for the single and double GDR 
states, respectively. As expected, the GDR peak is split in two while a broad 
plateau occurs in the DGDR region. The integrated cross sections are o"*^^^ = 
3.5 b and ct^gdr _ g g ^^ lower than those of Ponomarev et al. [13]. The 
difference can be ascribed to the fact that a coupled-channel method is used 
while in ref. [13] the cross sections are calculated in first- and second-order 



perturbation theory. In fact, if we use perturbation theory then the result for 
the (j*^°^ increases up to 4.2 b. The two approaches give similar results only 
for large impact parameters. Some difference comes also from the different 
-B(El) distribution. Figure 2 also shows the contributions associated with the 
0"*" (dashed line) and 2^ (dot-dashed line) component of the DGDR. The 1^ 
component does not appear in the figure since it is extremely small. After 
substraction of the long single GDR tail, the three peaks, expected from (10), 
are clearly visible. (If the convolution of the cross section were done with a 
r2 = 2ri, the three peaks are still visible though less evident). Therefore, an 
exclusive experiment in coincidence with 7-7 decay might conceivably give a 
direct signature of the excitation of the DGDR. 

The present results differ from those of Ponomarev et al. [13] who study the 
same reaction with the particle-phonon model, using second-order perturba- 
tion theory to calculate the cross section. In ref. [13] the DGDR cross section 
appears as a structureless peak while here it does not. The reason is that the 
calculated cross section to the single GDR in ref. [13] shows three peaks in the 
energy region 11-15 MeV. As a result, since the DGDR states are constructed 
as products of two single GDR states, one expects six peaks in the DGDR 
energy region which eventually smear out the cross section to the DGDR. 
In our case we have only two peaks, which are present in the experimental 
data, because we have fixed the parameters of our hamiltonian by fitting the 
photoabsorption cross section. 

The results of the calculated inelastic cross sections for the reactions ^"^^Pb 
-|- "^Sm at 0.5 GeV / A are shown in Fig. 3. The cross sections to each dis- 
crete state are shown as well as the ones obtained by the same smoothing 
procedure described above. In the transition from spherical to deformed one 
observes a continous evolution in the shape of the cross section resulting in 
the deformed case in a clear splitting in two peaks of the GDR and, corre- 
spondingly, three bumps in the DGDR energy region. This effect is due to the 
increase in separation between the two main components of the single GDR 
and the concentration of the small components in a more narrow energy range 
which in turn results from the increasing coupling of the GDR with quadrupole 
modes. 

In summary, energy and El properties of multiple giant dipole resonances have 
been studied in the context of the interacting boson model. The model has 
been applied to ^^^U and to transitional samarium isotopes for which Coulomb 
excitation cross sections have been calculated in the reaction with ^°^Pb at 
0.5 GeV / A. The example of the samarium isotopes shows how the excitation 
cross section is modified when going from spherical to deformed nuclei. The 
calculation shows also that exclusive experiments on a well-deformed nucleus 
like 238-[j could give a direct signature of the existence of the double giant 
dipole resonance. 
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Fig. 1. El excitation patterns up to (GDR)^ in (a) the spherical weak-couphng 
Umit and (b) the deformed strong-coupling limit, (a) The levels are labelled by J^ 
on top. The numbers between levels are the strengths -B(E1; (n, 0) J -^ (n + 1, 0) J') 
in units of C^ . (b) The levels are labelled by (A,^) underneath and 
by K"^ on top. The expressions between levels are the summed strengths 
Y.K'L> ^(El; (A,/x)KL -^ {X',n')K'L') in units of C^ 



Table 1 

Parameters of the hamiltonian and of the El transition operator. 

Isotope ep (MeV) ao (keV) 2^3 0^2 (keV) C (e fm) 
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Fig. 2. Coulomb excitation of the single and double GDR in ^ss^j jj^ ^^le 238 u ^q^ 
GeV/A) + 208p]-, pgaction. The dashed line is the contribution from the 0^ com- 
ponent of the two-phonon state while the 2^ one is represented by the dot-dashed 
line. The results corresponding to the 1+ component are too small to be seen. 
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Fig. 3. Coulomb excitation cross sections of single and double GDR in several samar- 
ium isotopes in the ^ospb (o.5 GeV/A) + Sm reaction. The bars correspond to the 
cross sections to the discrete states and the full line corresponds to the convolution 
of these cross sections by a lorentzian of width Fi = 2.5 MeV for the one-phonon 
states and r2 = V^Ti for the two-phonon states. 
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